Anisotropic ultrafast electron dynamics induced by high-field terahertz pulses in 

n-doped InGaAs 
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The anisotropic effective mass of electrons is directly measured using time-resolved THz- 
pump/THz-probe techniques in a ?i-doped InGaAs semiconductor thin film. A microscopic theory 
is used to attribute this anisotropy in the THz probe transmission to the nonparabolicity of the 
conduction band. Self-consistent light-matter coupling is shown to contribute significantly to the 
THz response. 
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PACS numbers: 73.50.Fq, 78.47 .jh 

Terahertz (THz) field induced transport effects can 
strongly influence the behavior of fast semiconductor 
components operating under extreme conditions [1, 2]. 
The internal dynamics of the charge carriers in such de- 
vices is, among other things, determined by the band 
structure of the solid. Consequently, a complete under- 
standing of the underlying physical processes requires 
a detailed knowledge of the band structure. In recent 
years, the development of extremely strong pulsed THz 
sources and ultrafast coherent detection methods [3-10] 
has made it possible to continuously alter the lattice mo- 
menta of the electrons in the solid and, thus, to explore 
their properties within the entire Brillouin zone. This can 
potentially yield new tools for characterizing the band 
structure and analyzing transport kinetics. 

Using various techniques of time-resolved THz nonlin- 
ear spectroscopy [4, 5, 10-26], effects like THz-induced 
intervalley scattering [15, 16] and ballistic transport of 
electrons across half of the Brillouin zone [21] have been 
investigated in semiconductors. Currently, also THz- 
pump/THz-probc (TPTP) techniques seem very attrac- 
tive in probing nonlinear carrier dynamics induced by 
intense few-cycle THz pulses in semiconductors. For in- 
stance, TPTP has been applied to induce impact ioniza- 
tion in InSb [14, 17] as well as intervalley scattering in 
doped GaAs, Si, and Ge [18]. 

For III-V semiconductor compounds, the conduction 
band energy £k is parabolic and symmetric only for low 
electronic momenta tik implying the same effective elec- 
tron mass in all directions, as illustrated in Fig. 1(a) 
at k = (0,0). As the electron distribution moves to- 
wards higher momenta at k = (K x ,0), however, e-^ be- 
comes strongly nonparabolic [27] producing strongly di- 
rectional, i.e. anisotropic, effective masses, as shown in 
Fig. 1(b) along the x- (black line) and y- (blue line) di- 
rections for a typical isotropic, nonparabolic band as a 



function of the corresponding electron energy. In other 
words, due to the nonparabolicity of £k, the two com- 
ponents of the effective mass arc equivalent only at the 
bottom of the band but they strongly deviate from each 
other at higher carrier momenta. 

Traditionally, the effective mass is measured using 
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FIG. 1. (color online) (a) The electric field of the THz pump 
pulse drives the electronic distribution (black circles) high in 
the band to an average momentum of k = (K x ,0), where 
the THz probe pulse senses the effective mass m x or m y in 
the x- or ^/-direction, respectively, (b) Effective mass nor- 
malized to the effective mass at the bottom of the band as 
function of the electron energy e Kx for co- (black line) and 
cross-linear (blue line) configurations. The effective mass is 
computed from Eqs. (2) and (3) using a nonparabolicity factor 
of a = 1.33 eV _1 . (c) Schematic of a polarization dependent 
THz-pump/THz-probe (TPTP) experiment. For the co-linear 
(CL) configuration, both THz pulses are polarized in the x- 
direction whereas the probe polarization is along the y-axis 
for cross-linear (XL) TPTP. 



2 




800 nm pump 



FIG. 2. (color online) Experimental setup, (a) Electric field 
profile of the THz pump beam emitted by the ZnTe optical 
rectification source. Inset: amplitude spectrum of the THz 
pulse, (b) Schematic of the experimental setup, (c) Electric 
field profile of the transmitted THz probe beam at various 
delay times between the main positive peaks of the THz pump 
and THz probe pulses. 



magnetic cyclotron resonance (CR) [28, 29]. In partic- 
ular, nonparabolicity was first observed in bulk InSb us- 
ing CR [28], and the same technique was used to study 
nonparabolicity in bulk Ino.53Gao.47As [29], showing an 
increase in the electron effective mass as a function of CR 
energy (up to 30 meV at 4.2 K). However, CR probes the 
average effective mass (often called cyclotron mass) for a 
given orbit in fc-space and for a certain level of energy, but 
is not capable of probing the above mentioned anisotropy 
of the effective mass at a specific point in fc-spacc. 

In this Letter, we present a new technique to directly 
monitor the anisotropic effective mass of the electrons. 
We implement a polarization dependent TPTP exper- 
iment as shown schematically in Fig. 1(c). Here, a 
strong THz pump pulse accelerates the electrons in the x- 
direction which are then probed by another weaker THz 
pulse polarized either in the x- (co-linear) or y- (cross- 
linear) direction. We show that the anisotropy in the 
electron masses yields distinctly different THz responses 
for the co- and cross-linear configurations because the 
measured THz probe signal is proportional to the in- 
verse effective mass in the probed direction. While the 
nonlinear THz response in doped semiconductors aris- 
ing from band nonparabolicities has already been re- 
ported [30, 31], our experimental results demonstrate the 
anisotropic nature of band nonparabolicity. 

In our experiments, a large aperture ZnTe optical rec- 
tification source [7] was used to generate high-power, 
few-cycle, THz pump pulses with energies of 0.8 \i J and 
0.1 — 3 THz bandwidth. Figure 2(a) shows an example 
of the temporal profile of the THz pulses produced by 
the ZnTe source, and the inset shows the corresponding 
amplitude spectrum. The TPTP setup is shown schemat- 



ically in Fig. 2(b). The sample used here was a 500 nm- 
thick n-type Ino.53Gao.47As epilayer (carrier concentra- 
tion of approximately 2 x 10 18 cm -3 ) on a 0.5-mm-thick 
semi- insulating InP substrate. A 10 x 10 x 0.5 mm 3 ZnTe 
crystal placed just after the first off-axis parabolic mirror 
and just before the sample was used to generate a THz 
probe beam that overlaps the THz pump beam at the 
focus on the sample. An additional ZnTe crystal 0.5 mm 
thick was used to detect the THz probe pulses transmit- 
ted through the sample by free-space electro-optic sam- 
pling. The spot size diameters on the sample for the 
THz pump beam and THz probe beam were 1.6 mm and 
2.5 mm, respectively (the probe beam path has a larger 
f- number than the THz pump beam). Both the ZnTe 
source crystal and the ZnTe detection crystal for the THz 
probe beam could be rotated to produce (and detect) 
probe polarization states either parallel or perpendicu- 
lar to the THz pump beam. At the sample position, the 
THz pump and THz probe peak electric fields were esti- 
mated to be about 200 kV/cm and 2 kV/cm, respectively. 
The THz probe itself is in the low-field regime and there- 
fore does not induce any nonlinear response in the sam- 
ple [15, 16]. We note that the non-collinear geometry of 
the TPTP experiment allowed the THz pump and probe 
beams transmitted through the sample to be geometri- 
cally separated. Cross-talk between the two THz beams 
was therefore avoided by simply placing a metallic beam 
block in the path of the transmitted THz pump beam 
after the second off-axis parabolic. In addition, lock-in 
detection of the transmitted THz probe pulse amplitude 
was synchronized to an optical chopper inserted in the 
THz probe source beam. As shown in Fig. 2(c), the am- 
plitude of the transmitted THz probe waveform increases 
when it overlaps with the THz pump pulse at zero rela- 
tive time delay, while the phase is relatively unaffected. 
This allows the transmission of the main positive peak 
of the THz probe pulse to be monitored as function of 
pump-probe delay time. We also note that all the ex- 
periments were performed under a dry nitrogen purge at 
room temperature. 

Figure 3(a) shows the normalized transmission of the 
main peak of the THz probe pulse as a function of pump- 
probe delay time, which is a common method for probing 
ultrafast carrier dynamics in semiconductors in optical- 
pump/THz-probe experiments [16, 32, 33]. The presence 
of the THz pump pulse results in an increase in trans- 
mission of the peak electric field of the THz probe pulse. 
The blue shaded area in Fig. 3(a) shows the transmission 
change for cross-linear polarization of pump and probe 
beams, while the black line shows the same measurement 
performed for the co-linear polarization configuration. In 
the latter fast, large amplitude oscillation is ob- 

served on top of a slower component similar to that shown 
as the area plot. As previously mentioned, we have ex- 
cluded the possibility that these fast oscillations are due 
to cross-talk between the pump and the probe at the 
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FIG. 3. (color online) Measured (a) and calculated (b) nor- 
malized peak transmission of the THz probe pulse as function 
of the pump-probe delay. The solid lines and the blue shaded 
areas show the results for co- and cross-linear TPTP configu- 
rations, respectively. The dashed line and the yellow shaded 
area in (b) represent corresponding calculations omitting the 
radiative backcoupling effects. 



detection level. In particular, when the probe beam is 
blocked, no residual signal from the pump beam is de- 
tected; both beams have to be present inside the sam- 
ple in order to observe the large amplitude oscillations 
for co-linear polarization shown in Fig. 3(a). Moreover, 
moving the sample to an off-focus position (or completely 
removing it from the THz beam) eliminated any modula- 
tion signal of the THz probe transmission. Finally, sim- 
ple interference cannot justify the observed anisotropy, 
due to the negligible contribution of the probe signal to 
the overall electric field amplitude. As discussed earlier 
for Fig. 2(c), the lack of any significant phase shift in 
the transmitted THz waveform implies that the observed 
signal is due to modulation of the peak amplitude of the 
THz probe pulse. To further corroborate this conclu- 
sion, no signal was observed as a function of pump-probe 
delay time when the THz detection point was set to a 
zero-crossing of the THz probe waveform. 

We next apply a microscopic theory to rigorously re- 
veal the microscopic origin of the detected anisotropy. 
The propagation of a THz field E(z,t) follows from the 
wave equation 
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where rib is the background refractive index, c is the 
speed of light, [i is the vacuum permeability, and J(i) 
is the excited current density. Since the thickness of the 
sample is much smaller than the THz field's wave length, 
the z-dependence of J is very accurately described via 
a 5-function. The temporal evolution of J follows then 
from J(t) = — j^j 5Z k (Vk£k) / k wriere /k defines the mi- 
croscopic electron distribution and V is the quantization 
volume of the semiconductor. For not too large fc-values, 
the conduction band energy dispersion is well-described 
by the relation [27] 
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where a is the nonparabolicity factor of the band and 
in* is the effective mass at the bottom of the band. The 
/c-dependent effective mass, i.e. 
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in Fig. 1(b) is computed using Eq. (2) with the standard 
InGaAs value a = 1.33 eV - [34]. To compute the dy- 
namics of the electron distribution, we use an equation- 
of- motion approach [35, 36] and finally obtain 
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The first term on the right hand side of Eq. (4) yields 
the acceleration of the electrons due to the THz field. 
This produces a time-dependent displacement of the elec- 
tronic distribution in k-space according to fk(t) = f^ t \ 

with 4zk(t) = — ^E, according to the classical accel- 
eration theorem [37]. The second term of Eq. (4) in- 
cludes Coulomb and phonon scattering of the accelerated 
electrons. This term leads effectively to a relaxation of 
the carriers and, thus, to a damping of the current den- 
sity. As we will discuss below, however, the experimen- 
tal results are almost not affected by these relaxation 
processes. Consequently, we can omit the corresponding 
term in the equation of motion. This reduces the numer- 
ical complexity and allows us to perform the calculations 
for a true three-dimensional system. 

Equation (1) shows that the THz field entering the 
equation for /k [Eq. (4)] itself depends on the electron 
dynamics via the current density. Consequently, Eqs. (1) 
and (4) must be solved self-consistently in order to ac- 
count for the backcoupling of the induced fields to the 
carrier dynamics. Physically, this backcoupling leads to 
a strong reflection of the THz field from the sample and 
to a radiative damping [36] of the induced current den- 
sity. For the conditions used in the experiment, we find 
effective radiative damping/relaxation times as short as 
50 fs implying a dominance of the radiative decay. Con- 
sequently, the results are only weakly affected by other 
scattering mechanisms. 
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We compute the transmission of a weak THz probe 
pulse in the presence of a strong pump field by numeri- 
cally solving Eqs. (1) and (4) for both the co- and cross- 
linear TPTP configurations. We start from a Fermi- 
Dirac distribution for /k at 300 K and an electron den- 
sity of 2 x 10 18 cm -3 . Similar to the experiment, we 
record the peak transmission of the probe pulse for var- 
ious pump-probe delay times. The results are presented 
in Fig. 3(b) and are in good qualitative agreement with 
the measured data [Fig. 3(a)]. We find that 3(t) pro- 
duces a high reflection of the THz probe pulse with only 
about 5% transmission in the absence of the pump pulse. 
At the same time, the magnitude of the current induced 
by the probe pulse is directly proportional to the in- 
verse effective mass. Therefore, the acceleration of the 
electrons into the nonparabolic regions can significantly 
modify J(t) and, thus, also the probe transmission de- 
tected. In particular, the solid line in Fig. 3(b) shows 
that the transmission is strongly enhanced when pump 
and probe are co-linearly polarized due to the increased 
effective mass. The cross-linear case (blue shaded area) 
produces a much smaller modification of the transmis- 
sion since the effective mass is changed only moderately 
[see Fig. 1(b)]. In both cases, the fast oscillations di- 
rectly reflect the transient pump-induced changes of the 
effective mass of the electrons. The theory analysis un- 
ambiguously assigns the difference of co- and cross-linear 
excitations to the mass anisotropy. The slow decay — 
only detected in the experimental signals — results from 
intervalley scattering [15, 16] which is not included in our 
model because it does not affect the detected anisotropy. 

The importance of the radiative backcoupling effects 
can be highlighted with an artificial calculation where 
we neglected these effects. Here, we simply replaced the 
full THz field in Eq. (4) by the external incident field. 
The corresponding results are shown in Fig. 3(b) as the 
dashed line (co-linear) and the yellow shaded area (cross 
linear). Evidently, such a simplification produces com- 
pletely opposite behavior than observed in the exper- 
iment. Although the basic findings of the experiment 
can easily be understood intuitively by the change of the 
effective masses, a complete understanding requires the 
inclusion of the complex self-consistent backcoupling ef- 
fects. 

In conclusion, we have used polarization-dependent 
THz-pump/THz-probe techniques to study the nonlin- 
ear THz response of electrons in n-doped InGaAs. As 
is unambiguously confirmed by our microscopic theory, 
the subpicosecond time resolution of our technique, cou- 
pled with the control of probe polarization, reveals the 
anisotropy in effective mass of hot electrons caused by 
the nonparabolicity of the conduction band. This new 
tool may open the way to directly mapping energy bands 
in semiconductors. Moreover, we have shown that the 
inclusion of backcoupling effects is crucial to understand 
the results. 
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